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We report a measurement of the top quark mass using events collected by the CDF II Detector 
from pp collisions at y/s = 1.96 TeV at the Fermilab Tevatron. We calculate a likelihood function 
for the top quark mass in events that are consistent with it —> bl~Vf,bl' + v'i decays. The likelihood 
is formed as the convolution of the leading-order matrix element and detector resolution functions. 

The joint likelihood is the product of likelihoods for each of 33 events collected in 340 pb _1 of 
integrated luminosity, yielding a top quark mass Mt = 165.2 ± 6.1(stat.) ± 3.4(syst.) GeV/c 2 . This 
first application of a matrix-element technique to tt —* M + uebl'~V(i decays gives the most precise 
single measurement of Mt in dilepton events. Combined with other CDF Run II measurements 
using dilepton events, we measure Mt = 167.9 ± 5.2(stat.) ± 3.7(syst.) GeV/c 2 . 

PACS numbers: 14.65.Ha, 13.85.Ni, 13.85.Qk, 12.15.Ff 


Precision measurements of the top quark mass, M t , 
place constraints on the masses of particles to which 
the top quark contributes radiative corrections, includ¬ 
ing the unobserved Higgs boson 1] and particles in ex¬ 
tensions to the standard model [21 At the Tevatron, 
top quarks are primarily produced in pairs. The dilep¬ 


ton channel, consisting of the decays it —> b£~Vebi ,+ 1 >' £ , 
has a small branching fraction but allows measurements 
which are less reliant on the calibration of the jet en¬ 
ergy scale, the dominant systematic uncertainty, than 
measurements in channels with hadronic W decays. A 
discrepancy from other channels could indicate contribu- 
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tions from new processes Q. 

The reconstruction of the top quark mass from dilepton 
events poses a particular challenge as the two neutrinos 
from W decays are undetected. Previous measurements 
in this channel Hi using Run I data have calculated the 
mass by making several kinematic assumptions and in¬ 
tegrating over the remaining unconstrained quantity. To 
extract maximum information from the small sample of 
dilepton events, we adapt a technique pioneered for the 
analysis of it —> Mvibqq' decays jfijilSlEGS]- This tech¬ 
nique uses the leading-order production cross-section and 
a parameterized description of the jet energy resolution. 
Making minimal kinematic assumptions and integrating 
over six unconstrained quantities, we obtain per-event 
likelihoods in top quark mass which can be directly mul¬ 
tiplied to obtain the joint likelihood from which M t is 
determined. 

This Letter reports a measurement with data collected 
by the CDF II detector between March 2002 and August 
2004, with an integrated luminosity of 340 pb -1 , yielding 
the most precise single measurement of M t in dilepton 
events, with statistical uncertainty approximately half 
that of Run I measurements 0 y|. We also report a 
combination of this measurement with other recent mea¬ 
surements using CDF Run II data in the dilepton chan¬ 
nel. 

The CDF II detector nn m is an azimuthally 
and forward-backward symmetric apparatus designed to 
study pp collisions at the Fermilab Tevatron. It con¬ 
sists of a magnetic spectrometer surrounded by calorime¬ 
ters and muon chambers. The charged particle tracking 
system is immersed in a 1.4 T magnetic field parallel 
to the p and p beams. Calorimeters segmented in 77 
and <f> surround the tracking system and measure the 
energies of interacting particles. The electromagnetic 
and hadronic calorimeters are lead-scintillator and iron- 
scintillator sampling devices. Drift chambers located out¬ 
side the central hadron calorimeters detect muons. 

The data are collected with lepton triggers that require 
events to have an electron or muon with px > 18 GeV/c. 
After offline reconstruction, we select events with (i) two 
leptons, each with px > 20 GeV/c, (ii) significant missing 
energy transverse to the beam direction (Jfix), and (iii) 
two jets, each with Ex > 15 GeV. The selection is de¬ 
fined as “DIL” in HU and was used to measure the cross 
section in the dilepton channel. These requirements yield 
33 events in the sample reported in this Letter. 

The probability density for it decays is expressed as 
P s (x\M t ), where M t is the top quark pole mass and 
x contains the measured lepton and jet momenta. We 
calculate P s (x\Mt ) using the theoretical description of 
the it production process expressed with respect to x, 
P s (x\M t ) = » where ^ is the differential 

cross section and a is the total cross section. 

To evaluate the probability density, we integrate the 


leading-order matrix element over quantities which are 
not directly measured by the detector, i.e., neutrino and 
quark energies. We assume that lepton momenta are 
perfectly measured, that quark angles are perfectly mea¬ 
sured by the corresponding jets, and that the two most 
energetic jets correspond to the b quarks from top quark 
decay. Quark energies, while not directly measured, are 
estimated from the observed energies of the correspond¬ 
ing jets. We define the transfer function W (p, j ) to be the 
probability of measuring jet energy j given quark energy 
p. We approximate W(p,j) as a sum of two Gaussians 
fitted to the predicted distribution of quark-jet energy 
difference from it events generated with Herwig and 
the CDF II detector simulation [J^. The expression for 
the probability density at a given mass for a specific event 
can be written as 

P s (x|M t ) = -^-y J d$\M t t{qi,Pi-,M t )\ 2 (1) 

x H W(pi,ji)fpDF(qi)fpDF(q2), 

jets 

where the integral is over the momenta of the initial 
and final state particles, q± and <72 are the incoming mo¬ 
menta, pi are the outgoing momenta, fpDFiqi) are the 
parton distribution functions and M t i(qi,Pi\M t ) is 
the tt production and decay matrix element as defined 
in Jjj,jl8j for the process qq —> it — > b(. + VfbP~Vp. While 
up to 15% of it pairs at the Tevatron are produced by 
gluon-gluon fusion ( gg — > it), this term can be excluded 
from the matrix element with negligible effect on the 
measurement. The term, l/a(M t ), in front of the in¬ 
tegral ensures that the normalization condition for the 
probability density, f dx P s (x|M t ) = 1, is satisfied. The 
effect of assumptions and approximations are measured 
in Monte Carlo experiments of fully realistic events, de¬ 
scribed below. Where needed, a correction is taken. 

We calculate the probability for the dominant back¬ 
ground processes, P bg (x), and form the generalized per- 
event probability density P(x|M t ) = P s (x\M t )p s (M t ) + 
P bgi (x)p bg 1 + P bg 2 {x)p bg2 • • •, where p s (M t ) and p bgi are 
determined from the expected fractions of signal and 
background events (see Table 0. The P bg are calculated 
in analogy to Equation 1, with the background matrix 
element evaluated numerically using algorithms adapted 
from the Alpgen jRfj generator. We calculate the prob¬ 
abilities for backgrounds arising from Z/ 7 * —> ee,pp 
plus associated jets, W+ > 3 jets where one jet is in¬ 
correctly identified as a lepton, and WW plus associated 
jets. Smaller backgrounds, comprising 11% of the ex¬ 
pected background, are not modeled. Application of the 
background probabilities reduces the expected statistical 
uncertainty by 15%. 

The posterior probability density in top quark mass 
is the product of a flat Bayesian prior and the joint 
likelihood, a product of the individual event likelihoods. 
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TABLE I: Expected numbers of signal and background events 
for a data sample of integrated luminosity of 340 pb _1 . Other 
backgrounds are negligible; all numbers have an additional 
correlated 6 % error from uncertainty in the sample luminosity. 


Source 

Events 

Expected tt (Mt = 178 GeV/c' 2 , 

a = 6.1 pb) 15.9 ±1.4 

Expected Background 

10.5 ±1.9 

Drell-Yan (Z/ 7 *) 

5.5 ± 1.3 

Misidentified Lepton 

3.5 ± 1.4 

Diboson (WW/WZ) 

1.6 ± 0.2 

Total Expected 

26.4 ±2.3 

Run II Data 

33 


The mass measurement (M t ) is the mean of the poste¬ 
rior probability, and the statistical uncertainty (A M t ) 
is the standard deviation. We calibrate the method us¬ 
ing Monte Carlo experiments of signal and background 
events. Signal events are generated with Herwig for 
top quark masses from 155 GeV/c 2 to 195 GeV/c 2 . 
Misidentified leptons are modeled using events from the 
data, while all other backgrounds are generated with 
Alpgen (Z/ 7 * —> ee, /r/r) or Pythia f2ol | (Z/ 7 * -*• 
tt, WW, WZ, ZZ) Monte Carlo. The numbers of sig¬ 
nal and background events in each pseudo-experiment 
are chosen according to Poisson distributions with mean 
values given in Table |T| The estimates for it yields 
take into account the mass dependence of the cross sec¬ 
tion [Si]] and acceptance. The response of the method 
for Monte Carlo experiments with both signal and back¬ 
ground is shown in Fig. ^ The response is consis¬ 
tent with a linear dependence on top quark mass, but 
has a slope that is less than unity due to the incom¬ 
plete modeling of background contributions. Corrections, 
M t -> 177.2 + (M t - 178.0)/0.84, A M t -+ AM t /0.84 are 
derived from this response and applied to both M t and 
A M t measured in data. 

Examining the width of the pull distributions in these 
Monte Carlo experiments, we find that the statistical un¬ 
certainty is underestimated by a factor of 1.51, indepen¬ 
dent of top quark mass. This results from simplifying 
assumptions described above, made to ensure the com¬ 
putational tractability of the integrals in Equation 1. The 
largest effects are from jets which originate from initial or 
final state radiation rather than &-quark hadronization, 
imperfect lepton momentum resolution, imperfect jet an¬ 
gle resolution, and unmodeled backgrounds. Correcting 
by this factor of 1.51, we estimate the mean statistical 
uncertainty to be 9.4 GeV/c 2 if M t = 178 GeV/c 2 @ 
or 7.8 GeV/c 2 if M t = 165 GeV/c 2 . 

We apply the method described in this Letter to the 
33 candidate events observed in the data. Including all 
corrections described above, we measure M t = 165.2 ± 
6.1(stat.) GeV/c 2 . Figure [ 2 ] shows the the joint prob¬ 
ability density, without systematic uncertainty, for the 
events in our data set. 



M t [GeV/c] 


FIG. 1: Mean measured Mt. in Monte Carlo experiments of 
signal and background events at varying top quark mass. The 
solid line is a linear fit to the points. 



FIG. 2: Joint posterior probability density as a function of 
the top quark mass for the 33 candidate events in data, after 
all corrections. Systematic uncertainties are not shown. 

The measured statistical uncertainty is consistent with 
the distribution of statistical uncertainties in Monte 
Carlo experiments where signal events with M t = 
165 GeV/c 2 are chosen according to a Poisson distri¬ 
bution with mean N t j = 21.7 events. This number of 
events corresponds to the cross section and acceptance 
at M t = 165 GeV/c 2 . Of these Monte Carlo experi¬ 
ments, 17% yielded a statistical uncertainty less than 6.1 
GeV/c 2 . 

There are several sources of systematic uncertainty 
in our measurement which are summarized in Table UH 
The largest is due to the uncertainty in the jet energy 
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TABLE II: Summary of systematic uncertainties. 


Source 

A M t (GeV/c 2 ) 

Jet energy scale 

2.6 

Limited background statistics 

1.2 

PDFs 

1.0 

Generator 

0.8 

Background modeling 

0.8 

FSR modeling 

0.7 

ISR modeling 

0.5 

Response correction 

0.4 

Sample composition uncertainty 

0.3 

Total 

3.4 


scale ‘23j, which we estimate at 2.6 GeV/c 2 by vary¬ 
ing the scale within its uncertainty. An uncertainty of 
1.2 GeV/c 2 comes from the limited number of back¬ 
ground events available for Monte Carlo experiments. 
The contribution from uncertainties in the PDFs are es¬ 
timated by using different PDF sets (CTEQ5L 0 vs. 
MRST72 |H), different values of Aqcd, varying the 
eigenvectors of the CTEQ6M 0 set, and varying the 
initial state contributions of gg and qq ; the total corre¬ 
sponding uncertainty added in quadrature is 1.0 GeV/c 2 . 
Dependence on the Monte Carlo generator is estimated 
as the difference in the extracted top quark mass from 
Pythia events and Herwig events; this amounts to 
0.8 GeV/c 2 . We observe no difference in the extracted 
top quark mass in events from a leading-order and a next- 
to-leading order generator 0 ■ We estimate the un¬ 
certainty coming from modeling of the two largest sources 
of background, Z/y* and events with a misindentified 
lepton, to be 0.8 GeV/c 2 . The uncertainty due to imper¬ 
fect modeling of initial state (ISR) and final state (FSR) 
QCD radiation is estimated by varying the amount of 
ISR and FSR in simulated events |2fj and is measured 
to be 0.7 GeV/c 2 for FSR and 0.5 GeV/c 2 for ISR. The 
uncertainty in the mass due to uncertainties in the re¬ 
sponse correction shown in Fig. □ is 0.4 GeV/c 2 . The 
contribution from uncertainties in background composi¬ 
tion is estimated by varying the background estimates 
from Table [I] within their uncertainties and amounts to 
0.3 GeV/c 2 . Adding all of these contributions together 
in quadrature yields a total systematic uncertainty of 3.4 
GeV/c 2 . 

Applications of other dilepton techniques [28] using the 
same CDF Run II data yield values consistent with this 
result. The four results are combined using a standard 
method [29’]. We determine statistical correlations be¬ 
tween the measurements using Monte Carlo experiments, 
and assume systematic uncertainties to be 100% corre¬ 
lated except the few that are method-specific which are 
assumed to be uncorrelated. Table mu gives the four 
dilepton measurements, their statistical correlations, and 
their weight in the combination. The NWA method uses 
looser event selection and thus has a smaller statisti¬ 


cal correlation. Correlations significantly less than unity 
suggest that each method extracts unique information. 


TABLE III: Measurements of Mt in the dilepton channel with 
statistical and systematic errors, their statistical correlations 
and the weight of each measurement in the combined result. 


Method 

Result (GeV/c 2 ) 

Correlation Matrix 

Weight 

This Letter 

165.2+^] ±3.4 

1 



0.47 

NWA [28] 

170.7±6 5 ±4-6 

0.12 

1 


0.36 

KIN [28] 

169.51^2 ± 4.0 

0.40 

0.14 

1 

0.18 

PHI [28] 

169.7±9.o ± 4.0 

0.43 

0.25 

0.35 1 

0.00 


Combining the four CDF Run II dilepton measure¬ 
ments we obtain 

M t = 167.9 ± 5.2(stat.) ± 3.7(syst.) GeV/c 2 . 

This combined result is consistent with the current world 
average (30| and the sin gle m ost precise measurement in 
the lepton+jets channel |3ll| . 
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